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Abstract. Two models of the circumstellar disc around supergiant B[e] stars are discussed: an equatorial wind 
model produced by wind bi-stability, and a Keplerian viscous disc model. Both models are successful in providing 
a site for dust formation once they have cooled sufficiently. However, the optical-infra-red continuum is calculated 
and it is found that both models have significant trouble in accounting for observations. In particular the optical- 
near-IR emission is accounted for, but the dust emission is underestimated by at least an order of magnitude. 
Variations in the structure of the models (the temperature variation with radius, the density structure and the 
dust opacity) are investigated to assess how (in)appropriate the standard models are for supergiant B[e] star discs. 
Changing the temperature structure, and making simple dust opacity changes within the disc has little effect on 
the resultant continuum emission. By altering the density structure of the discs, the continuum may be accounted 
for by both models: the equatorial wind model requires a very flat density profile which is impossible to explain 
with any accelerating wind, and the viscous disc model’s density structure is required to fall off less steeply with 
radius than would have been expected, although this may be explained from consideration of viscous processes in 
the disc. It is recognised that both theoretical interpretations have difficulties and unsolved problems. 
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1. Introduction 

Observations of supergiant B[e] stars (denoted sgB[e] stars 
hereafter - see Lamers et al. 1998 for a classification sys¬ 
tem) show several features: infra-red (IR) continuum ex¬ 
cesses and far-IR dust emission (e.g. see Zickgraf 1992); 
broad (~ 10 3 km s _1 ) UV resonance lines; in the opti¬ 
cal there are strong Balmer lines (Ha equivalent widths 
of > 10 2 A), and narrow (10s of km s _1 ) permitted and 
forbidden line emission. These stars have typical stel¬ 
lar parameters of M* ~ 30-70M Q , i?* ~ 30-100 -Rq , 
T e s ~ 15,000-30,000K (Zickgraf et al., 1986). 

Zickgraf et al. (1985) proposed that the observations 
could be explained with a multicomponent wind struc¬ 
ture: a fast line-driven wind (appropriate for a hot star) 
over the poles and a dense disc of outflowing gas in the 
equatorial plane. This geometry has received some cor¬ 
roboration from polarimetry (Zickgraf & Schulte-Ladbeck 
1989, Oudmaijer et al. 1998). 

Whilst there is a good physical explanation of the fast 
polar wind, (Castor, Abbott & Klein, 1975) the formation 
mechanism (and indeed the structure) of the dense equa¬ 
torial flow is still largely unknown. An excellent attempt at 
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explaining the equatorial flow has been made by Lamers & 
Pauldrach (1991 - rotationally-induced wind bi-stability), 
which has been recently advanced by Pelupessy, Lamers 
& Vink (2000) who suggest that wind compression due to 
its rotation (Bjorkman & Cassinelli 1993) may also play a 
part in the generation of the disc (although wind compres¬ 
sion was found by Owocki, Cranmer & Gayley 1996 and 
also by Petrenz & Puls 2000 to be strongly inhibited in 
rotating winds). Zickgraf et al. (1996) fit a wind-like flow 
with a “beta” type velocity law to observations of the UV 
metallic lines of R50 and Hen S22 (interpreted as edge-on 
stars) and found acceptable fits with wind terminal veloc¬ 
ities of 60-80 km s _1 . 

Oudmaijer et al. (1998) present a hydrodynamical 
model of the flow around HD87643 in which a line-driven 
wind is driven partly by luminosity from a Keplerian re¬ 
processing disc. As well as producing a fast polar wind, 
their model produces a region several scale-heights above 
the equatorial plane which has intermediate velocities 
(~ 100 km s _1 ), and hence they point to this region being 
responsible for the permitted line emission. However, they 
did not discuss the disc formation, being more concerned 
with the intermediate and fast radiatively driven flows. 

Given the current uncertainties in the sgB[e] stars’ disc 
origin, a model for an equatorial wind produced by the 
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bi-stability mechanism, and also a model for a viscous 
Keplerian disc (similar to that currently gaining accep¬ 
tance in the classical-Be star community) are directly com¬ 
pared. The aim of this study is to assess if either model 
can be ruled out, or (at least) if one model can be strongly 
favoured over the other as the more likely disc structure 
of sgB[e] stars. As a first step, the optical IR spectral 
energy distribution is calculated and compared to obser¬ 
vations (the line emission is returned to in the discussion). 
This is not a trivial task as a successful model must also 
provide a site for the formation of dust. In §2 the density 
and velocities of both the equatorial wind model and the 
viscous disc model are discussed, along with a heuristic 
argument for the temperature profile. The possibility of 
the flow forming dust is considered in §3. In §4 the contin¬ 
uum emission is calculated and its variation with relevant 
parameters presented. A discussion is given in §5 and con¬ 
clusions in § 6 . 

2. The models 

The main concern of this study is to examine if either 
model can account for the continuum observations of 
sgB[e] stars. Most previous work has examined aspects 
of wind models (for example Bjorkman, 1998, presents 
results of a wind compression model) whilst this is the 
first time a viscous disc model has been applied to sgB[e] 
stars. During the course of the study, values for the 
stellar parameters, fast polar wind etc. are taken from 
the Large Magellanic Cloud star R126: M* « 40M 0 , 
i?* « 7 Oi ?0 , T « 22, 500K (Zickgraf et al. 1985), the 
fast polar wind parameters vary widely in the litera¬ 
ture: a mass-loss rate of M w = 1 O _ 6 - 1 O _ 5 M 0 yr _1 , and 
terminal velocity of Voo « 1800 km s -1 was calculated 
by Zickgraf et al. (1985), whilst Bjorkman (1998) gives 
M w = 4.6 x lO” 5 M 0 yr _1 and Voo = 650 km s -1 from 
Kudritzki et al .'s (1989) “cooking recipe”. The adopted 
nomenclature used in this paper is that mass-loss rates 
will be quoted as a mass-flux multiplied by 47 t i.e. the rate 
is quoted as if the mass-loss was present over the whole 
star. This should be particularly noted when considering 
the mass flux through the disc. 


terised by a latitudinal varying mass flux rate M(8) and 
wind speed v r (0). At the jump latitude 6 , M(9) increases 
and v r decreases producing a large increase in density 
p(6) oc M(Q)/r 2 v r . Within this high density, low veloc¬ 
ity region (representing the “disc”), the radial velocity is 
v rt d and density is pd■ The model requires that the star 
is rotating at a substantial fraction (> 50%) of its criti¬ 
cal speed (where centrifugal forces balance gravity). The 
azimuthal velocity of the wind in the disc decreases 
with radius from the star due to angular momentum con¬ 
servation and so is not very important dynamically in the 
disc. The disc thickness is determined by the value of 6 
where the bi-stability jump occurs, and so depends on the 
rotation rate and how close the effective temperature is to 
the jump temperature. 

Within the disc, the density and radial and azimuthal 
velocities follow the form 


Pd 


Vr,d 


^4>,d 


M d 


47T r 2 v r ,d 
w(r)voo,d 



(1) 

( 2 ) 

(3) 


where v^o is the rotational speed at the stellar surface 
and again note the disc mass-loss rate is quoted as if this 
mass loss applied over the whole star - the actual mass- 
loss rate through the disc is M d multiplied by the solid 
angle of the disc divided by 47 t. The velocity law is quoted 
in the form above with w(r) = (1 — Ro/r where i?o = 
(1 — and Wi is the fraction of the terminal velocity 

that the flow has at the star. It is expected that the radial 
velocity at the stellar surface (corresponding to toj) will be 
less than the sound speed. With terminal velocities Vqo ~ 
100 km s _1 (Zickgraf et al. 1985) and the sound speed ~ 
10-20 km s -1 for B stars this corresponds to vj. t « 0.1. 
Zickgraf (1992), estimates the disc opening angle to be 
10-20° from the relative frequency of sgB[e] stars with 
edge-on and polar characteristics, and this value is used 
hereafter. 


2.1.2. Viscous disc model 


2.1. Density and velocity structure 

2.1.1. Equatorial wind model 

The popular model for the equatorial density enhance¬ 
ment of sgB[e] circumstellar media is one based on 
radiatively-driven winds. Support for this geometry came 
in 1991 when the bi-stability model was presented, in 
which the stellar wind parameters (mass-loss rate and ve¬ 
locity) “jump” over a band centered on the equator of 
the star as the driving species and lines change at tem¬ 
peratures in the region of 25,000K due to the changing 
conditions at the photosphere from a gravity-darkened 
fast rotating star (Laniers & Pauldrach 1991, Vink et al., 
1999, and Pelupessy et al., 2000) The model is charac¬ 


Tlie model of viscous outflowing discs has been discussed 
by Lee et al. (1991) and Okazaki (2001) and applied to 
classical Be star observations by Porter (1999). In this 
model, the angular momentum distribution within the disc 
is determined by viscosity, in a similar way to the more 
familiar accretion disc. However, in this case the angu¬ 
lar momentum source for viscous discs is the central star 
(which acts as an angular-momentum sink for accretion 
discs). In common with the bi-stability mechanism, this 
model also requires that the star is a fast rotator - indeed 
it insists that the equatorial region is rotating at Keplerian 
speeds. The weak link in the model for classical Be stars 
(which will also feature in this application to sgB [e] stars) 
is that the mechanism to inject angular momentum in to 
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the disc at its inner regions is still largely unknown - often 
quoted candidates are pulsation and magnetic fields. 

The solution of the equations of mass continuity and 
angular momentum conservation in the radial and az¬ 
imuthal co-ordinates produce a disc which has surface- 
density, radial and azimuthal velocity are well represented 
in the inner ~ 10 2 J?* by power laws in radius, (Okazaki 
2001 ). 

In particular the discs are close to Keplerian v$ oc 
r^ 1 / 2 , the radial velocity increases linearly with radius 
ty oc r (and is subsonic until large radii), and the surface- 
density E (the density integrated through the disc) de¬ 
creases as E oc r~ 2 . The disc flares - the scale height H 
of the disc increases with radius as H = c s r 15 /s/UMl 
where G is the gravitational constant, M* is the stellar 
mass and c s is the sound speed. The viscous disc density 
and velocities are 


Pd = P0,d ^ 

^R* j 


( r 

Vr,d = VrO,d 

V R* 

( 

r \ 

V<t>,d = Vqo 1 

R * ) 


- 3.5 


exp 


- 0.5 



(4) 

(5) 

( 6 ) 


where po,d, VrO,d and v^o are the density, radial velocity, 
and azimuthal velocity at the star-disc boundary (r = 
R * ) respectively, and z is the height above the equatorial 
plane. Also, v r0 ,d <, 10 ~ 2 c s and v^o « \JGM * /i?* In the 
outer parts of the disc, it becomes transonic, and angular 
momentum conserving (changes to ~ r^ 1 ). 


2.2. Temperature profile 

The form of temperature profile for a disc is difficult to 
determine without a full calculation. However a simple 
argument can be applied to both equatorial winds and 
viscous discs which can provide an estimate of the radial 
variation of temperature. It is assumed that the disc will 
have a constant temperature with radius until a cooling 
radius r c , and beyond this the disc will cool. In the limit 
that no radiation of any wavelength impinges on the disc, 
it will cool adiabatically such that T oc r -4 / 3 . However, it 
is most likely that the disc will act as a reprocessing disc: 
it will absorb radiation from the star, and re-emit at its 
local equilibrium temperature Td- Hence the total energy 
emitted by the disc (oc T t f) is proportional to the flux re- 
cieved, which in turn is proportional to r~ 3 (r~ 2 from the 
dilution of the radiation field, and an extra factor of r _1 
from the cosine term produced by the geometric reduction 
of the projected disc area). Hence these reprocessing disc 
have T oc r -3 or oc r -3 / 4 . This argument is valid for 
both equatorial winds, and viscous discs. 

The cooling will change the metallic species responsible 
for the absorption and scattering of the equatorial wind 
and hence may lead to a more complicated discussion of 
the wind driving. This has already been implicitly taken 


into account by using the slowly accelerating beta-velocity 
law fitted to the UV lines from Zickgraf et al. (1996). 

Cooling will effect the viscous disc: once the tempera¬ 
ture in the viscous disc varies, then the density and veloc¬ 
ity profiles are also affected. This is due to the torques 
changing as the viscosity v = ac s H alters (a is the 
Shakura & Sunyaev 1973 viscosity parameter). This is a 
result of the sound speed c s and also the disc scale height 
H oc c s changing. To ensure mass and angular momentum 
conservation, the density, radial velocity and scale height 
change to 


p oc r" 10 / 8 
v r oc r 1 / 4 
H ex r 9/8 


(7) 

( 8 ) 
(9) 


i.e. their radial dependence becomes less severe, whilst the 
azimuthal velocity remains Keplerian (Okazaki 2001, pri¬ 
vate communication). Therefore at radii where the viscous 
disc is cooling (r > r c ) eqn.7 9 replace eqn.4-6. 


3. Dust formation 


Any model which attempts to explain the properties of 
sgB[e] star circumstellar matter must also provide a site 
for dust formation. To produce dust, two criteria need to 
be fulfilled: the temperature of the gas must be lower than 
the sublimation temperature of the dust (around 1500K 
depending on the chemical composition of the dust), and 
secondly, the number density of the species involved in the 
formation of the dust needs to be above a critical value. 

It is expected that sgB[e] stars’ dust will be silicate 
based, as carbon is depleted due to the CNO process¬ 
ing cycle in massive stars, and indeed sgB[e] stars do 
show features due to silicates in their spectra (e.g. Voors, 
1999). Examination of the 8-13/rm spectra of R126 pre¬ 
sented in Roche et al. (1993) shows that there are no 
strong emission features, which suggests silicate dust is 
not present. However, suppression of the silicate feature 
may be achieved by having large grains (ad > 10/im), or 
by not having a stratified disc structure (leading to any 
optically thin “atmosphere” region) i.e. the emitting re¬ 
gion is completely optically thick (see Dullemond et al, 
2001 and Meeus et al. 2002). 

Gail & Sedlmyar (1988) showed that a critical num¬ 
ber density can be derived by comparing the timescale for 
the chemical reactions responsible for dust formation r c h 
to the timescale for a parcel of gas to change in density 
and temperature r eX p (i.e. an expansion timescale). Grain 
growth will occur if r c h < r eX p. Gail & Sedlmyar find that 
this occurs for any size of grain if the number density of 
the growth species is larger than a minimum density 


^min 


V r 

10 -16 ri> re i 


( 10 ) 


where u re i is the relative velocity between the growing 
grain cluster and a molecule (set to the thermal veloc¬ 
ity here). It is assumed that dust formation occurs if 
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Fig. 1. Comparison of number density of carbon (solid 
line) with n m in from eq.10 (dotted line): top panel, the 
equatorial wind model described in the text, lower panel, 
the viscous disc model. Dust may form in the region where 
the number density is larger than ?i m ; n only if the flow is 
cool enough. The arrows indicate the radius at which the 
temperature has fallen to 1500K. 

ep(r)/ (fj,m p ) > n m i n (where the ratio of relevant species 
to total number density is e and p « 1.25 -1 is the mean 
molecular weight for neutral gas). This criterion may 
be compared to differing models of the disc component 
around sgB[e] stars. 

The value of e is taken to be e = 10 -5 correspond¬ 
ing to 1 /3 of Solar silicon abundance appropriate for the 
LMC metallicity. As an example calculation of n n i m , the 
sound speed (approximating to w re i) has been calculated 
assuming that for both the equatorial wind and the vis¬ 
cous disc models, cooling takes place beyond r c = 1.5i?* , 
and that the temperature profile is T oc (r/r c ) -3 / 4 there¬ 
after. Finally the temperatures at r = i?* are taken to be 
0.8T e g and 0.5T e g' for the equatorial wind and disc models 
respectively. The results of this calculation are shown in 
fig-1. 

3.1. Equatorial wind model 

The equatorial wind is assumed to have the same param¬ 
eters as described in §2.1.1: (/ 3 = 4 and w t = 0.1) and a 
mass-loss rate of M ( j = 1O -5 M 0 yr _1 (this rate has been 
chosen as it produces a dense enough disc to acccount for 
the near-IR emission of the model star R126 (see §5). Fig.l 
shows the number density of this wind model as well as 
the critical number density required to form dust. It can 
be seen that dust can only form for r < 501?* . The arrow 
on fig.l shows the radius where the temperature falls to 


1500K - this occurs at a radius of 41i?* . Therefore there 
exists a region from ~ 40-50??* where the gas is dense 
enough and cool enough for dust formation. 

3.2. Viscous disc model 

Fig.l also shows the number density within a vis¬ 
cous disc with po = 10~ n g cnr -3 (again chosen to 
reproduce the near-IR emission, see §5) and v r o = 
0.3aT 4 ^R* /M* km s _1 (where a = 0.1 is the viscosity 
parameter, T 4 is the disc temperature in 10 4 K and ??* and 
M* are in Solar units - see eq. 6 of Porter, 1998). Again 
the cooling radius is set to r c = 1.51?* . For r <; 12501?* 
the number density in the disc is high enough to enable 
dust to form. The radius where the temperature falls to 
1500K is 23f?* , well within the region where dust may 
form. 

3.3. Dust formation conclusions 

Both the viscous disc model and the wind model have a 
region where their number density is larger than the crit¬ 
ical value for dust formation, implying that both models 
can form dust, with the proviso that the temperature must 
also be below the sublimation temperature Tdust- The ma¬ 
jor difference between the models is for n m in, which is sig¬ 
nificantly lower for the viscous disc model as its radial 
velocity is far lower than for the equatorial wind model. 
Both models are successful in providing a site for dust 
formation, and so both are still good candidates for the 
actual discs od sgB[e] stars. 


4. Continuum emission 

The continuum emission is now calculated from the den¬ 
sity fields in the models and compared to observations 
of R126. The model includes a fast polar wind (with a 
fixed mass-loss rate of lO _5 M0yr _1 ), the disc models 
described in the last section, and dust (where it exists) 
in the disc. Although the fast polar wind is included, its 
contribution to the final flux is negligible, and is only in¬ 
cluded for completeness. In order to calculate the excess 
emission in the optical - IR, the prescription of Waters 
(1986) is used. The star is assumed to be pole-on. 

The viscous disc model requires that the density po,d, 
the dust formation temperature Td us t and the cooling ra¬ 
dius r c are specified to calculate the emission. The equato¬ 
rial wind model requires a mass-loss rate a terminal 
velocity Voo^d, the value of /? for the velocity law and a 
flow opening angle, as well as the dust temperature and 
cooling radius: f3 is fixed at f3 = 4, the terminal velocity 
Voo = 70 km s-\ and the opening angle to 20° hereafter. 

The number density is calculated where the tempera¬ 
ture falls below T^us t and compared to the minimum num¬ 
ber density for dust formation from eq.10. If dust is found 
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to form, the optical depth through the disc at that radii 
is 

Tdust = / ptvdustdz (11) 

J disc 

where p is the gas density, and tc^ust the opacity of the 
dust. The opacity used is an approximate fit to that cal¬ 
culated for a distribution of grain sizes (e.g. see Wood et 
al. 2002). A fit to Wood et aids data for Solar metal- 
licity in their fig.2 from 0.3-100/im produces Kdust = 
120(A/0.1/xm ) -0 ' 6 cm 2 g~ 1 . The opacity used here is 1/3 
of this, to take into account the lower metallicity of the 
LMC. Therefore the dust opacity used here is 

“ d “ = 40 (ou^) cm2g "‘- (12) 

This simple expression for dust opacity is used so that 
the reason for the failure of the standard viscous disc and 
equatorial wind models - see below - is clear. In order that 
the dust is readily observed in the spectrum, it is expected 
that the optical depth of the dust will be Tdust > 1 - 

The emitted flux is then integrated over the disc area, 
which, for face-on discs, produces 

1 /*^disc 

F v = J TrB u [T(r)](l - e~ Tduat )2m'dr (13) 

where i?disc is the outer edge of the disc, and D is the 
distance to the star and B v [T(r)} is the Planck function. 



Fig. 2. Dust optical depths within the discs at a wave¬ 
length of 5/xrn for the two “standard” models. The solid 
line is the log of the optical depth (referring to the left 
hand scale), the dashed line is the log of the disc scale 
height H in stellar radii (again referring to the left hand 
scale), and the dotted line is the log of the disc density 
in g cm -3 (the right hand scale). The top panel is for the 
equatorial wind model, the bottom panel for the viscous 
disc model. 


5. Results 

The UBVRIJHKLMN photometry for R126 is taken from 
Zickgraf et aids (1985) compilation, and is supplemented 
by 12/zm and 25pm IRAS points from Schwering (1989). 
The reddening to the object was taken to be E(B- 
V) = 0.25 from Zickgraf et al. (1985). The calculated emis¬ 
sion from the model has been reddened with the expres¬ 
sion given in Howarth (1983) and the model flux in the 
V band has been normalised to the R126 observations for 
all models. The underlying photospheric continuum is a 
Kurucz 22000K, \ogg = 3.0 ATLAS9 model. 

The proceedure used to fit the photometry involves 
a search over the mass-loss rate M,j. dust formation 
temperature and cooling radius for the equatorial wind 
model. The viscous disc model has the density normal¬ 
ization poi dust formation temperature and cooling ra¬ 
dius varied to find the best fit. The results are sum¬ 
marised in Table 1: the various parameters for the mod¬ 
els are given, along with their evaluation of fit. Best- 
fit values to the optical-near-IR data produces values of 
Mr/ = 1.1 x 10 - 5 M©yr -1 and po = 10~ lo ' 75 g cm -3 for 
both models respectively. The equivalent mass loss rate 
for the viscous disc model is 7 x lO~'M 0 yr _1 : although 
it appears counterintuitive that this is less than the po¬ 
lar (fast wind) mass loss rate, it is reconciled in that the 
mechanism for mass ejection into the viscous disc is not 


the same as over the poles (indeed the cause for mass injec¬ 
tion in to the disc it is still unknown). Hence the equatorial 
mass loss rate need not bear any relation to the polar mass 
loss rate. 

There is a significant problem with the dust emission 
for both models in that the optical depths where dust 
forms is not large enough to generate the observed emis¬ 
sion. To illustrate this, fig.2 shows the dust optical depth 
Tdust (and density and disc scale-height) for both models 
as a function of radius. It can be clearly seen that the 
optical depth is significantly less than unity at the dust 
formation radius, and decreases with radius. Therefore it 
might be expected that the continuum emission will not 
be able to reproduce the obserevational data. Fig.3 dis¬ 
plays the full spectrum for both models. As can be seen 
the optical to near-IR flux is reproduced, although long- 
ward of ~ 2 pm the dust emission is too low by about an 
order of magnitude. 

How may this be remedied? Or are the two proposed 
models to be rejected? To produce the observed emission, 
the optical depth Td us t needs to be increased by one to 
two orders of magnitude. For these models this requires 
a significant increase in density. A global increase in the 
density (corresponding to increases in M,j and po) would 
produce too much emission in the optical-near-IR range, 
and so may be ruled out. Increasing the opening angle 
for the equatorial wind model (or decreasing the terminal 





Fig. 3. Continuum emission from the star (R126) and the 
two best-fit models. The total model flux (solid line) con¬ 
sists of Kurucz model atmosphere (long dashed line), free- 
free and free-bound emission (dotted line) and dust emis¬ 
sion (short-dashed line) for the “standard” models. The 
model has been normalised to the R126 data of Zickgraf 
et al. (1985) in the V band. The left panel is for the equa¬ 
torial wind model, and the right panel is for the viscous 
disc model. Model parameters are given in Table 1. 

velocity) will also produce the same effect of increasing 
the emission at all wavelengths, and not simply the dust 
emission. To investigate how dependent the models are to 
the input assumptions, several changes are made in turn to 
the standard models: the temperature power law indices 
are made free parameters, the dust opacity is set to that 
appropriate to Solar metallicity, and finally, the density 
power law indices are made free parameters. 

5.1. The temperature profile 

It is possible that the temperature profile that has been 
used is not correct, although the disc must cool in order for 
it to form dust. To assess whether the observations may 
be fitted using a different temperature profile, beyond the 
cooling radius, the power law Td oc r~ m has been varied 
using different values of m, from 0.4 to 1.5 and the best 
fit models calculated. The best-fit parameters are shown 
in Table 1. The temperature-power law index of the best 
models is to = 0.6 for the equatorial wind and to = 0.5 
for the viscous disc. 

In the case of the equatorial wind the best-fit dust tem¬ 
perature is identical to the standard wind model (1400K), 
and the cooling radius is slightly smaller. However, the 
fit is again characterised by a significant underestimate of 
the dust emission indicating that variation of the tempera¬ 
ture power law produces a negligible effect on the resultant 
spectrum. The viscous disc produces a slightly cooler dust 
formation temperature (1500K) than in the standard case 


Fig. 4. As fig.3 except that the temperature power law in¬ 
dex is a free parameter. The left panel is for the equatorial 
wind model with T oc r -0 , and the right panel is for the 
viscous disc model T oc r~ 0 ’ 5 . Further model parameters 
are given in Table 1. Comparison with fig.3 illustraltes 
that the effect of allowing the temperature power law to 
change is negligible. 

(1600K). However, as in the wind case, the fit is barely 
improved with significant underestimate of the dust emis¬ 
sion. 

Both of these models are presented in fig.4. It appears 
that the failure of the standard models to reproduce the 
observations is not dominated by the temperature profile 
in either model. 

5.2. The dust opacity 

Whilst so far the opacity from Wood et al. (2002) has been 
used, corrected to take account of the lower metallicity of 
the LMC, here the normalization for Solar metallicity is 
used (i.e. the LMC star R126 is assumed to have Solar 
metallicity). This illustrates the effect of increasing the 
opacity, whilst keeping the form of the density field the 
same as in the standard cases. 

The dust opacity used here is ftdust = 
120(A/0.1/zto) - 0 - 6 cm 2 g -1 , and the best-fitting models 
are shown in fig.5 (with parameters in Table 1 : models 
5 and 6). Formally, from the % 2 values in Table 1, it can 
be seen that increasing the dust opacity decreases the 
residuals. However, the best fitting models (fig.5) are still 
woefully inadequate in accounting for the dust emission. 

Further increase in the dust opacity should provide 
a high enough optical depth to account for observations. 
However, this would require an exceptionally large metal¬ 
licity for this star which is physically very difficult to jus¬ 
tify at this stage of evolution of the star, and therefore 
this possibility is disregarded in favour of the dynamical 
arguments given later. 
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Table 1 . The parameters for the best-fit models. The first two models are the standard models with m fixed at 0.75. 
The following two models have the temperature power law index m as a free parameter in the search and with the 
inner flow parameters as above. Models 5 & 6 have an enhanced dust opacity appropriate for Solar metallicity. The 
final two models have a power law in density as stated - see text. 


Model 

m 

r c /R* 

Tdust (K) 

Kdust (0.1pm) 

1. Wind 

0.75 

2.3 

1400 

40 

2. Viscous disc 

0.75 

2.0 

1600 

40 

3. Wind 

0.60 

1.1 

1400 

40 

4. Viscous disc 

0.50 

2.1 

1500 

40 

5. Wind 

0.75 

2.0 

1300 

120 

6. Viscous disc 

0.75 

2.9 

1400 

120 

7. Wind 

0.75 

3.5 

700 

40 

8. Viscous disc 

0.75 

6.5 

800 

40 


5.3. The density profile 

The standard density profiles for the wind model (eq.l) is 
now changed to a pure power law form: p/po = ( r/R * )~ n 
(where po is a normalization density at r = A* ). The 
sole difference between the two models now is (i) that the 
equatorial wind model has a disc with a constant opening 
angle (20°), whereas the viscous disc flares, and (ii) that 
the density falls exponentially with height away from the 
equatorial plane for the viscous disc. Only the temperature 
power law index is fixed (with m = 0.75), and the best-fit 
model search is conducted over dust temperature Xd us t, 
density power law index n, density normalization po, and 
cooling radius r c . 

The best-fitting equatorial wind model is shown in 
fig.6, and provides an excellent description of the contin¬ 
uum flux of R126. The parameters of this fit are n = 1.7 
and po = 10 _11 ' 9 g cm -3 , cooling radius of r c = 3.5A* 
and dust temperature of T<j us t = 700K. A literal interpre¬ 
tation of this best fit solution would predict that the radial 
velocity actually decreases with radius (from comparison 
with eq.l), which indicates that the simple wind model is 
an erroneous description of the actual disc around R126. 

The best fiting viscous disc model is also shown in fig.6 
and also provides a good fit to the data. The density of 
this model is po = 10 _11 °g cm -3 (note this is significantly 
higher than the best fitting equatorial wind model), and 
the density power law index is n = 2.7, the cooling radius 
is r c = 6.5 and the dust formation temperature is Tiust = 
800K. 

For both the equatorial wind, and the viscous disc the 
y 2 (n, po , r c , Tdust) values have been inspected for differ¬ 
ent models around this minimum, in order to estimated 
the variation in the parameters. A value of y 2 is obtained 
which is twice the best fitting model if parameter uncer¬ 
tainties are Sn ~ ±0.1, Spo ~ ±0.1 ST^ust ~ ±100A', and 
Sr c ~ ±3. Clearly the least constrained parameter is the 
cooling radius. 

6. Discussion 

It appears from the results that both the “standard” equa¬ 
torial wind model and the viscous disc model have trou- 


density y 2 


M w 

= 1.1 x 10 -5 M q yr _1 

407 

P = 

10- lo ' 75 (r/7?.* )- i5 g cm’ 3 

403 

M w 

= 1.1 x lO^Moyr" 1 

359 

P = 

Kr la7 (r/A* )- 3 ' 5 g cm" 3 

397 

M w 

= 1.1 x lO^Moyr -1 

316 

P = 

KT 10 - 7 (r/A, )- 3 ' 5 g cm" 3 

313 

P = 

10 _11 ' 9 (r/A* ) _1 ' 7 g cm" 3 

64 

P = 

KT 110 (r/A* ) _2 ' 7 g cm' 3 

61 



log(A/ gm) log(A/ fan) 

Fig. 5. As fig.3 except that the dust opacity normal¬ 
ization is appropriate for Solar metallicity: ftdust = 
120(A/0.1/im)~°- 6 g cm -3 . Further model parameters are 
given in Table 1. The left panel is for the equatorial wind 
model and the right panel is for the viscous disc model. 


ble in generating the continuum emission for sgB[e] stars. 
To do so, the density field has had to be changed some¬ 
what from initial theoretical expectations. The deviations 
from the standard density structures of the models may 
indeed highlight some aspects of the two models which 
have not been considered here. The equatorial wind may 
be clumped at large radii, hence producing parts of the 
disc which have a large dust optical depth. However, the 
total emission may not necessarily increase due to the sub¬ 
sequent reduction in effective emitting area (i.e. if the gas 
is clumped into dense regions, the “filling factor” of that 
dense gas degreases). The low power law index for the best 
fitting viscous disc model (n = 2.7 compared to n = 3.5 for 
the standard model) might indicate an extra dependence 
of the viscosity (and perhaps the a parameterization) on 
the surface density, or the disc temperature which has not 
been investigated here. 
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0.5 0 0.5 1 It®.5 0 0.5 1 1.5 

log(\/ /im) log(X/ H 


Fig. 6. As fig.3 except that the density power law is a free 
parameter. The left panel is for the equatorial wind model 
with p = 10 -11 ’ 9 (r/i?* ) _1 ' 7 g cm -3 , and the right panel is 
for the viscous disc model p = 10 _ll o (r/i?* ) _2 7 g cm -3 . 
Other model parameters are given in Table 1. 


Assuming that either of these models can account for 
all the continuum emission, can the two models both pro¬ 
duce the observed heavy element line emission from sgB [e] 
stars? The equatorial wind model was introduced to do ex¬ 
actly that (e.g. see Zickgraf et al. 1996), and the line emis¬ 
sion comes from the equatorial wind flow itself. Therefore, 
this model would seem to provide a complete description 
of the disc around a sgB[e] star. However, there is a dif¬ 
ficult issue for the equatorial wind model: the electron 
scattering optical depth becomes larger than unity at rel¬ 
atively small radii. Whilst high continuum optical depths 
are commonplace for (approximately) spherically symmet¬ 
ric Wolf-Rayet winds, in this asymmetric case it may lead 
to problems in driving the wind. Photons cannot physi¬ 
cally penetrate into the disc which leads to little or no 
radiative driving for the gas, and hence the density profile 
may not actually allow the gas to be driven! This might 
cast doubt on whether the equatorial wind model is cor¬ 
rect in principle. However, this may also highlight that the 
wind may not have a velocity law typified by the “beta” 
parameterization. 

The production of the heavy element line emission 
from the viscous disc model is slightly more complicated. 
The dominant velocity component in the disc is a ro¬ 
tational one, and therefore optically thin emission lines 
should exhibit double-peaked profiles (similar to those 
observed in classical Be stars). However, the model in 
Oudmaijer et al. (1998), provides a separate site for line 
emission with a dominant radial component of velocity: 
a disc wind. Here, the emission from the reprocessing 
disc helps to drive a flow from the upper and lower sur¬ 
faces of the disc to relatively modest terminal velocities 
~ 100 km s _1 . The resultant line emission may have a 


double-peaked disc component, but will have a large ra¬ 
dial component imposed on it too. Currently, no attempt 
has been made to synthesize line profiles from the work 
presented by Oudmaijer et al. An issue to resolve here is 
whether there is enough flow generated in the disc-wind 
to account for the lines observed in sgB[e] stars. 

A weakness for the viscous disc model is that the 
source of angular momentum required to supply the disc 
is still unidentified (although candidates exist: pulsation 
and/or magnetic fields). Until a coherent theory of angular 
momentum transfer from the star to the disc is produced, 
then the viscous disc model will only remain a promising 
candidate. Note that if the viscous disc model is the cor¬ 
rect structure, then it is possible that the way in which 
classical Be stars supply angular momentum to their discs 
is the same for sgB[e] stars. In this case then the sgB[e] 
stars will be invaluable as they have a different envelope 
structure than classical Be stars, but it might be likely 
that both supply angular momentum in the same fashion. 

From the results above, it is difficult to make a strong 
case which will prefer one model over the other, or even 
that either model is appropriate. Clearly, the discs of 
sgB[e] stars are not the simple radiation driven flows, nor 
the simple viscous discs that have been previously sug¬ 
gested. 

7. Conclusion 

Both of the models that have been presented in this study 
are apparently incapable of reproducing the observational 
results in their simplest form. Whilst both meet the 
requirements for the formation of dust they are unable to 
yield the necessary broad-band fluxes. 

The first model is the enhanced equatorial flow from 
Lamers & Pauldrach’s bi-stability theory (1991); the other 
is the viscous disc model. The densities are alike, but the 
velocity fields are disparate. Both competing paradigms 
expect that the central supergiant star is rotating at a 
significant fraction of its break-up speed. 

By relaxing several assumptions (regarding the tem¬ 
perature profile, the dust opacity, and the density pro¬ 
file), the applicability of the models have been tested. It is 
found that the models are unable to reproduce the obser¬ 
vations unless the disc density field is significantly altered, 
although then a good representation of the observations 
can be obtained. 

There are possibilities for both models to produce the 
observed permitted and forbidden line emission: the equa¬ 
torial wind model produces the emission itself, and the 
viscous disc model produces a disc-wind from it’s sur¬ 
faces which may produce the lines (although this remains 
untested). 

In both the equatorial wind model and the viscous disc 
model, there is still work to be undertaken to understand 
the form of the density profile necessary to account for 
observations. Clearly both models in their simplest form 
are inadequate as a description of the discs of sgB[e] stars. 
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